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Abstract 

We show that the gluon distribution function for very large nuclei may be computed 
for small transverse momentum as correlation functions of an ultraviolet finite two 
dimensional Euclidean field theory. This computation is valid to all orders in the 
density of partons per unit area, but to lowest order in a s . The gluon distribution 
function is proportional to 1/x, and the effect of the finite density of partons is to 
modify the dependence on transverse momentum for small transverse momentum. 
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1 Introduction 



In a previous paper, we argued that in a limited range of transverse momentum, 
for small values of Bjorken x, quark and gluon distributions functions for very 
large nuclei might be evaluated as the solution of a weakly coupled many body 
theory Q. This result relied heavily on the technology of light cone quantization. [Q] 
- H Specifically, when a parameter 

corresponding to the density of charge squared fluctuations per unit area, is fj? » 

Aq CD , then the strong coupling is a s (y?) « 1. When the Bjorken x is x << A^ 1 / 3 , 

it is then valid to replace the valence quarks by delta functions of charge along the 

light cone. For a s /u 2 << q 2 << fj, 2 , and in this range of x, we computed the gluon 

distribution function to lowest order in weak coupling to be 

1 dN _ a s y 2 (N 2 -l) 1 
irR 2 dxd 2 qt vr 2 xq 2 

The gluon distribution function per unit area was precisely the Weizsacker-Williams 

distribution function for gluons scaled by the density of charge squared fluctuations 

per unit area. The physical picture corresponding to the above formula is that the 

Weizsacker-Williams distribution is generated by random fluctuation in the charge 

per unit area, and is similar in spirit although different in origin than pictures used 

to describe nucleus-nucleus scattering. ||-[]|] 

The approximation of small x guaranties that the central region gluons see a 

source of valence quarks which are of a much smaller size than a typical gluon 

wavelength as measured in a frame co-moving with the gluon distribution at that 

value of Bjorken x. We are therefore in the deeply screened region. In this kinematic 

region, Lipatov enhancements of the gluon distribution function are expected to 

modify the Bjorken x dependence of the distribution function and take 1/x — ► 

l/ x l+Ca s where Q 

is some constant ||. If such a small x enhancement actually 
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occurs, then the weak coupling expansion which is allowed will only be formal since 
a s {n 2 ) ln(l/x) will not be small. It is not yet clear whether this enhancement 
actually takes place in the deeply screened small x region we are interested in. If 
it does, although the coupling constant is weak, one will have to find a way of 
systematically including the effects of this enhancement. 

We will not address the small x enhancement in this paper. We shall instead 
turn to another aspect of the problem-which is, computing the distribution functions 
in the small qt region. We will here work to lowest (formal?) order in a s but to all 
orders in ct^/i 2 . We will show that the correlation function which gives the gluon 
distribution function can be expressed as a two dimensional Euclidean correlation 
function of an ultraviolet finite field theory. We will find that summing to all 
orders in a 2 /x 2 modifies the qt distribution function. However, to all orders in this 
expansion, the gluon distribution function is proportional to 1/x. 

Before deriving these results, we shall first briefly review the results of our 
earlier paper. We recall that in the small x region, the valence quarks are Lorentz 
contracted to a smaller distance scale than that of the wavelength of the co-moving 
gluon. Therefore the valence quarks may be treated as being approximately delta 
functions along the light cone. We can ignore valence quark recoil, so long as the 
gluons being emitted are soft and so long as the coupling is weak. In this limit, the 
valence quarks may be treated as sources of charge. 

The valence quarks may also be taken to be uniformly distributed in trans- 
verse space for sufficiently large nuclei. In this case, the scale of variation of the 
nuclear valence quark distribution along the transverse direction can be made large 
compared to the typical hadronic distance scales. 

The problem therefore is to compute the distribution function of gluons in the 
presence of static sources of color charge localized along the light cone and uniform 
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in transverse space. The external current due to the source may be represented as 



It was shown in our first paper that the nucleus can be broken up into regions 
of transverse spatial extent such that the number of valence quarks in each region 
is large. This allows us to treat the sources of charge as classical. We also showed 
that the dominant contribution to the ground state wavefunction came from states 
which had large fluctuations around zero color charge, but where the fluctuations 
were small compared to the total amount of charge in each transverse spatial region. 
In this limit, the fluctuations are gaussian. 

The A 1 / 3 dependence of the gluon distribution function in this screened region 
follows from the above arguments. While the color charge is screened so that the 
average color charge is of order y/~N, where N is the total color charge in each spatial 
region, the coherence of the color field makes the gluon density of order TV. The 
number of valence quarks per unit transverse area goes as A 1 / 3 . The gluon density 
therefore goes as yl 1//3 -up to corrections due to the logarithmic dependence on A of 
the coupling constant. 

The problem is therefore a simple one: If we want to compute a ground state 
correlation function, we can do it by the path integral 



that is, we just integrate the path integral for fixed charge around a Gaussian 
fluctuating charge at each point in space. We must therefore compute correlation 
functions in a stochastic background field. 

The approximation that we may treat the source as classical is only true in the 
limit where the spatial regions we are looking at have a large number of quarks in 
them. In our first paper, we argued that this requires that 



x t )S(x ) . 



(3) 




(4) 



Qt «V 2 



(5) 
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We will assume that this is also the case in this paper. 

In the previous paper, we evaluated the ground state gluon distribution function 
in the perturbative region where both a s and a 2 /j, 2 were treated as small parameters. 
This latter condition forced a 2 pL 2 « q 2 . In this paper, we will relax this condition 
and set up the computation of the gluon distribution function in the soft region. 

2 The Classical Problem 

We first turn to the problem of computing the solution of the classical problem for 
the gluon field in the presence of a source which is a delta function along the light 
cone. The equation of motion is 

D li F^ = 9 r, (6) 

where J is the classical light cone source. We will work in light cone gauge where 
A_ = -A+ = 0. 

The current J M only has components along the + direction and is proportional 
to a delta function of x~ . There exists a solution of the equations of motion for 
this problem, where the longitudinal component A+, which is not zero by a gauge, 
vanishes by the equations of motion 

A+ = -A~ = 0. (7) 

The only non-zero components of the field strength therefore are the transverse 
components which we require to be of the form 

Mx) = 9(x-) ai (x t ) , (8) 

where 6(x~) = +1 for x~ > and 9{x~) = for x~ < 0. This function is 
nonvanishing only when x~ > 0, which since x~ = (t — x)/y/2 is equivalent to 
x < t. This is what we expect for a classical field generated by a source traveling 
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close to the speed of light with x = t. For x > t, the source has not yet arrived, 
and for x < t the source should produce a field. 
Using the definition of F l+ in terms of A 1 

F i+ = 5{x~)ai. (9) 

If we further require that 

F ij = 0, (10) 

(where i and j are transverse components), we see that we have a solution of the 
equations of motion so long as 

V-a = gp(x t ). (11) 

Here p is the surface charge density associated with the current J. There is no 
dependence on x~ because we have factored out the delta function. The dependence 
on x + goes away because the extended current conservation law, 

d+Q a + if abc A b + Q c = 0, (12) 

is simplified by the solution of the field equations in Eq. (7) to read 

8 + Q a = 0. (13) 

Hence, 

Q a (x + ,x t )=p a (x t ). (14) 

The condition that i™ = is precisely the condition that the field a is a gauge 
transform of the vacuum field configuration for a two dimensional gauge theory. 
The requirement that V • a = gp, is a condition that fixes the gauge. For such a 
field configuration, the light cone Hamiltonian 

P~ = 0. (15) 



This is the precise analog of what we found in our previous paper for the Weizsacker- 
Williams distribution around an electron. 

The field configuration which is a gauge transform of the vacuum field configu- 
ration for a two dimensional field theory may be written as 

T-Oi = — U- ViC/t. (16) 
9 i 

We have not been able to construct explicit solutions for the above equation for 
arbitrary dependence of the surface charge density on xt- 

Note that the field configuration which solves this problem only has a nontrivial 
dependence on x t . The dependence on x~ , as shown in Eq. (8) is only through a 
step function, and upon Fourier transforming gives only a factor of l/k + . As we 
will see in the next section, the distribution functions associated with this field 
therefore, to all orders in a 2 s y?, have the general form 

>?/*■>■ 

In our previous paper, we showed explicitly that in the weak coupling limit we 
obtained the simple result 

H(y) = l/y, (18) 

for 

l«y«l/a$. (19) 

We see here that this result is true only in the weak coupling limit where a^fj? « 
kf. The range is larger than that used in our previous work. Note also that this 
result is parallel to that of finite temperature field theory where there is a non- 
perturbative length scale, the magnetic screening length, which is A ~ l/a s T where 
T is the temperature. 
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3 Computing Correlation Functions 

To compute correlation functions associated with our classical solutions, we must 
integrate over all p. This is equivalent to integrating over the transverse field with 
the constraint that the field must be a pure gauge, that is, 



We can perform the integration over the transverse field a as integrations over 
unitary matrices in the standard way in which one goes to a latticization of a gauge 
theory. The integration measure becomes 



We see that the measure for this theory is that for a two dimensional Euclidean 
field theory. The spatial variables are those of the two dimensional transverse space 
of the original theory. This theory is ultraviolet finite because of the fact the La- 
grangean is fourth order in derivatives. The expansion parameter for the theory is 
is a 2 s p 2 /k 2 . 

The above analysis ignores Fadeev-Popov determinants which are generated in 
transforming from the integration over p to the integration over a and finally to that 
over U. To properly define the two dimensional theory for purposes of Monte-Carlo 
simulation, one may have to investigate these determinants. For our purpose, which 
is to study the scaling behavior of expectation values, we will not take into account 
these determinants. The determinant for both integral representations in terms of 
a or in terms of U is easy to compute, and is the same for both representations. 
It is the determinant of V ■ D where D is the covariant derivative associated with 
the field a. This leads in two dimensions to an ultraviolet finite modification of the 
above measure. 

As discussed in the previous section and more explicitly in our earlier paper, 
when the relevant momentum scale is k 2 » a^n 2 , the theory is in the weak coupling 




(20) 




(21) 
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region and may be evaluated perturbatively . For smaller values of k 2 , we are in the 
strong coupling phase of the theory. In this phase of the theory, we expect that there 
should be no long range order. Correlation functions of xt should die exponentially 
at large distances, or alternatively the Fourier transform of correlation functions 
should go like a constant for small momentum. We will see that this guarantees the 
finiteness of the gluon distribution functions for small momentum. 

It should be easy to compute correlation functions for the above action using 
lattice Monte Carlo methods. The theory is two dimensional which should make 
possible the use of large grids. The theory is Euclidean so that all quantities of 
physical interest are computable as Euclidean correlation functions. The theory is 
ultraviolet finite, so that there should be no problems extrapolating to the contin- 
uum limit. 

Finally, the correlation function for the computation of the transverse momen- 
tum dependence of the structure function must be determined. This may be simply 
evaluated to be 

D{k t ) = —H{k t ) (22) 

where D is the propagator for the two dimensional theory 

(2ir) 2 5 2 (k t - q t )5 tj D{k t ) =< a t a 3 > (23) 

The expression for the propagator can be easily rexpressed in terms of the link 
variables U. As claimed, the small kt behavior of the gluon distribution function 
is related to the asymptotics in coordinate space for the propagator of the two 
dimensional theory. Large distances correspond to strong coupling, which in turn 
corresponds to a lack of correlations signalled by an exponential fall off. We therefore 
expect that at small k t , the function H(kt) will be finite; the gluon distribution 
function will be non-singular. 
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4 Summary and Conclusions 

We have seen that the infrared behavior of the gluon distribution function for small 
kt is determined by solving a two dimensional field theory. The relevant fields 
are gauge transforms of two dimensional vacuum configurations. The light cone 
Hamiltonian vanishes for such configurations. The field theory for the correlation 
functions of interest is finite and involves integrating over gauge transforms of the 
field with a specified weight function. 

It is clear that the gluon propagator and light quark propagator in the presence 
of such a background field configuration should be quite simple. The solutions of 
the small fluctuation equations are just gauge transforms of free field solutions. To 
construct the propagator one must join solutions across the discontinuity in x~ 
generated by the source of charge. The study of these propagators will be the 
subject of later analyses. 

It is also clear that the infrared behavior of these correlation functions is com- 
putable as a lattice Monte Carlo simulation. Realistically it appears that the desired 
accuracy might be obtained. 
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